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Abstract. We solve the problem of the temporal evolution of one of two-modes 
embedded in a same dissipative environment and investigate the role of the losses 
after the preparation of a coherent state in only one of the two modes. Based 
on current cavity QED technology, we present a calculation of the fidelity of 
a superposition of coherent states engineered in a bimodal high-Q cavity. Our 
calculation demonstrates that the engineered superposition retains coherence for 
large times when the mean photon number of the prepared mode is on the order 
of unity. 

Problems involving interaction of atoms and two modes of a electromagnetic field 
trapped into a high-Q cavity have increasingly received attention [1-14]. From the 
theoretical point of view, bimodal cavities were considered in a number of papers. For 
example, in the phenomena of a degenerate and nondegenerate two-mode squeezing 
for a generalized Jaynnes-Cummings model with a two-level atom was considered; in 
the authors show how to prepare W states, GHZ states, and two-qutrit entangled 
states using the multi-atom two-mode entanglement; in 0E1E1 schemes for teleporting 
entangled state of zero- and -one photon state from a bimodal cavity to another were 
proposed; in the authors show how to build bilinear and quadratic Hamiltonians, 

thus opening the possibility to implement up- and down-conversion operations in two 
mode cavity QED. Recently, experiments with coherent control of atomic collision and 
controlled entanglement in bimodal cavities have been reported [§] , attracting even 
more attention to this topic. 

It is to be noted that all the above mentioned papers do not take into account the 
effect of the environment, which plays an important role mainly when the decoherence 
time has to be considered. Damping of two modes of a cavity QED has been considered 
in Refs. [121 1131 114) . all of them using the Liouville approach. In the behavior of a 
two-level atom interacting with two modes of a light field in a cavity was investigated, 
and the effects of cavity damping were treated numerically for the special case of 
fields initially prepared in two mode squeezed vacuum state. In ref. |13| the authors 
studied the evolution of a two-level atom coupled to two modes of a dissipative cavity, 
and in ref. ^3] the authors have investigated cross decay rates and robust (against 
dissipation) states in a system composed by two cavity modes in the presence of the 
same reservoir. 



Superposition of coherent states prepared in one mode. 



2 



In this paper we review the problem of two modes in a lossy cavity QED. Our 
model differs from that of J3| since we take into account bosonic modes in a dissipative 
environment. Also, different from ^1 our model includes, beside two non- 
interacting bosonic modes subject to the same reservoir, a two-level atom interacting 
off-resonantly with one of the two modes of the bimodal cavity. Our Hamiltonian 
model is 



H = H + Hi (1) 



where 



and 
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H = ^2 + E ^kblh + —^-cr z + ha\aiXicr z , (2) 



H i = E E n ( A **°5 6 * + A >4) • (3) 
k j=i 

Here we consider a two-level atom composed by ground \g) and excited |e) state, 
°~z — l e ) ( e l — Iff) (fli ! Qj an d dj are, respectively, the creation and annihilation operator 
for the ]th cavity mode of frequency u)j, whereas b\. and bu are the analogous operators 
for the kth reservoir oscillator mode, whose corresponding frequency and coupling with 
the mode j = 1, 2 write tOk and Xjk, respectively. The atom- field coupling parameter is 

2 

Xi = where g is the Rabi frequency and 5i = (wi — wq) is the detuning between the 
field frequency u>i and the atomic frequency luq. Note that the dispersive interaction 
occurs with only mode 1. This is important for preparing superposition of coherent 
states PH]- Using the completeness relation for both atom and fields as given by 
coherent states, the Schrodinger state vector associated with Hamiltonian Eq.Q can 
be written as 

Mt)) = \g)\Mt)) + \e)\Mt)) (4) 

where |^(t)) = / ^ / £sa J |^!&} A t \a u a 2 , {f3 k }), I = g,e. The complex 

quantities ai,ot2,{0k} stand for the eigenvalues of a,j and bk, respectively, and 
Ai (ai, a2, , t) are the expansion coefficients for \<pe(t)) in the basis of the coherent 
state products |ai, a%, The assumption of |{/3&}) for the reservoir follows 

from a remarkable property first proved by Mollow and Glauber |15j . i.e., that at 
zero absolute temperature the reservoir receives coherently the excitation lost by a 
coherent state. The extension to finite temperature is considered in Ref.|16). where 
a phenomenological-operator approach to dissipation in cavity QED was proposed. 
Essentially, this extension consists in treating the reservoir plus the surrounding as 
a pure state provided that we include the degrees of freedom of whatever remain in 
the surrounding, in such way that we obtain a mixed (thermal) state after tracing out 
the degrees of freedom of whatever was included. Because the orthogonality of the 
atomic states and Eq. Ill II) , we can obtain the uncoupled time-dependent Schroedinger 
equations 

ihj t \Mt)) = He\Mt)) , (5) 

where Hi = Hi + Hj , and 

H m = h(wi± x) a \ a i + h\oia\ai + ^ huj k b\bk- (6) 
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Note that the problem has been reduced to that of dissipation of two modes of the 
cavity field whose frequency w\ of mode 1 has been shifted by — x (x) when interacting 
with the ground (excited) state of the atom. To solve the problem of the evolution 
of an arbitrary initial state prepared in one of the two modes we pursue a different 
approach from that of Refs. |12U13llT^| . Instead of following the evolution of the two 
modes at once, we will follow rightly to the mode of interest, assumed as mode 1. A 
convenient way to perform this can be done by the method of the reduced density 
operator ^2]. For this purpose we first calculate the characteristic function G in the 
normal order, then the Glauber-Sudarshan P representation, and finally the reduced 
density operator p\{t) for mode 1. The G function for mode 1 in the Heisenberg 
picture reads 

G(v,V*,t) = Trm {p UR (0) exp na\{t) exp H7*oi(*)]}, (7) 

where pi2n(0) is the density operator for the whole system composed by modes 1,2 
and the reservoir at the instant t = 0, and Tr 2 R indicates the trace on mode 2 and 
reservoir. The Glauber-Sudarshan P representation is given by the two-dimensional 
Fourier transform of the characteristic function G: 

P(t,e,t) = J Gfo, + »7*0^r. (8) 

while the reduced density operator for mode 1 of the cavity field state is given by the 
diagonal representation of P (£, £*, f): 

pi(t) = fp(t,t,t)\Z)(Z\<Pt- (9) 

To apply this method, we will need to solve the following Heisenberg equations for 
modes dj corresponding to the Hamiltonian Eq.JHJ) in order to substitute them in 
Eq.0: 

dj = - iujiOj - X* k b k (10) 
fc 
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b k = - iu k bi - «2J Ajfctti, (11) 
i=i 

where ojji — (u>i ± x) if j = 1 and the minus (plus) signal depends on the initial state 
g (e) of the atom, or u>jt — u)% if j = 2. Eas. HlUlll|) can be solved, for example, by 
Laplace transform. In doing so, we multiply both sides of Eas. (|10lllfl by exp(— st), 
integrate from to oo and evaluate the inverse Laplace transform under the Wigner- 
Weisskopf approximation ^21 

s + tuJk 2 

where Aw^-, jjj denote the shift in the energy and the decay rate of the j-mode, 
respectively, and Aujji, jjf (j j') the corresponding cross shift in the energy and 
the cross decay rates |T3] for modes 1,2. The results for modes 1 and 2 are 

ai(t)=«n(t)oi(0)+ui 2 (t)a 2 (0)+5^1?i it (t)6 fc (0) ) (12) 

fc 

a 2 {t)=U22(t)a2(0)+U2i{t)a 1 {0) + J2^k{t)b k (0), (13) 
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un(t) = exp 



+ cosh 



Ui2{t) = - exp 



and 



(A + B) 



(B-A) 



(B- A) Z +A\CY 



: sinh 



(B - A) 2 + 4 |C|" 



(B — A) 2 + 4 | C | 2 



(A + B) 



2C 



_^{B-A) 2 +A\C\ 



■ sinh 



\b-ay + a\c\' 



(14) 



(15) 



1 ? lfc (t) = i(A|E 2 C-AiE 1 -AJ 1 B) 
exp 



x < 



+ 



A+B-^/{B-A)^+i\C\ 2 
2 



exp 



^+B+ A /(B-A) 2 +4|C|") ^ 



( ( b + a) + v(b-^) 2+ 4| ct _ Wfe)v / (B _ A)2+4 |cr 

exp(— iujt) 



t (B+A) + ^(B-Ay>+A\C\ 2 _ ^ -(B+A) + ^/{B-Af+A\C\ 2 
\ 2 luJ k)\ 2 



MJfc) 



(16) 



with 



A = i (lui ± x + Awi) + 7n/2 (17) 
B = i(w 2 ±X + Aw 2 ) +712/2 (18) 
C = - tAwu - 712/2. (19) 

Functions U2i(£), U22W and $2k(t) can be obtained from wi2(i), wn(i) and dik{t), 
respectively, simply replacing C by C* , A by -B, 1 by 2 and vice-versa. 

It is now straightforward to obtain pi(£) for an arbitrary initial state if we note 
that any initial state J2 n m Cnm \m) 1 \n) 2 can be lwritten in its most general form as 



Pi2fl(0) = ^2 d 2 a l d 2 a 2 d 2 (3 1 d 2 f3 2 C nmjk \ (a 2 ||/3i) 22 (/3 2 | 

n,m,j,k 

x J {d 2 J k }P({l k })\{l k }) R (b k }\ (20) 

where modes 1 and 2 were expanded in coherent states and the reservoir is in a 

thermal state J^Pk \ n k) r ( n k\ characterized by an average photon number in the kth 

k 

mode (rife) = 1/ [exp(fiuk/kT — 1)] according to pt = , _[" fc L + i , and was expanded in 
the diagonal representation with P ({7/0}) = Yl x (n fc ) ex P ("fej) ■ ^he characteristic 
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function for mode 1 corresponding to the initial state of Eq. (|2U[1 reads 
G(rj,r]*,t)= ^ J d 2 aid 2 a 2 d 2 (3id 2 f3 2 C nrnj k 2 (A2I /3i) 2 i ("2! oc x ) l 

n,m : j,k 

x cxp [-rfu\x{t)a. i + Wiii^a^} ex P [-f?*ui2(*)/3oi + rju^t)^] 

x I] cx p[-M 2 |<Mi)| 2 (n*)] , (21) 

where we have written /(0) = /o- The Glauber-Sudarshan P t) representation 

for mode 1 as given by Eq. iJHJ is the Fourier transform of Eq. (|21fl : 

(&| Pi) 2 



P(£,e,t) = J d 2 a 1 d 2 a 2 d 2 /3 1 d 2 p 2 C nmjk (f3a) 1 



n,m,j,k 



nD(t) 



( \j - mi(t)a i - u 12 (t)/3 w ] [C - ut^afo - u* 12 (t)P* 20 } \ , . 
V D(t) ) K 1 

where D(t) = \$ik{t)\ 2 (nk)- The reduced density operator for mode 1 can now be 
obtained using Eq.©. 

In the following, let us consider the important case of engineering two coherent 
states inside a single bimodal high— Q cavity with one of the two modes supporting a 
coherent state and the other mode being prepared in a superposition of coherent state 
(SCS) in the presence of a reservoir cooled to K, i.e., when the initial state is given 

by 

K (HO))! + \-<*(0))i) H°)>2 m) R ■ (23) 
This example is the analogous of the "Schrodinger cat state" studied in Ref.^H] m 
the context of unimodal cavity QED, in the microwave domain, to unravel the role 
of the environment in the transition from the quantum to classical dynamics. For 
calculating the corresponding evolution of the SCS of Ea. l23|) . it is enough to consider 
the evolution of the initial state Pi 2 r(0) as given by Eq.J^OJ), and left a\, a 2 and /3i, (5 2 
assuming the values ±a to compose the complete density operator. Using Eas. (|21l22|l 
and taking the limit D(t) — > for zero temperature, we obtain for the evolved mode 
1 in Eq.JU 

^ ,^ 2 (fho\ (3io) 2 1 (ago] Quo)i 

1 (-itii(t)a!2o + U2i(t)Pox\ un(t)a w + u 12 {t)(3 02 ) l 
x \u 11 {t)a 1Q +u 1 2{t)Poi) ll {-u n {t)a2Q+uvi{t)P()2\- (24) 
Now, taking into account Eq. l|24|l . when considering the initial state as given by 
Eq.lO, Pi(t) reads 

pi(t) = Af e {\[u\i{t) + ui 2 (t)] a w } 1 ([uix(t) + ui 2 (t)] aio\ 
+ + ui 2 (t)} a w ) 1 ([-un(t) + ui 2 (t)] a w \ 

+ Z(t)[|[un(t)+ui2(*)]a 1 o) 1 ([-uii(t)+u 12 (t)]oio|+/i.c.]} (25) 
where h.c. means Hermitian conjugate and 

2 Mi 2 (t)Mn(t) u* 12 (t)un(t)' 



Z(t)= exp<j-2|a | 2 l-|« u (t)| 



(26) 



is the term responsible by decoherence. To evaluate the errors in the SCS prepared 
in one of the two modes caused by the environment and due to the presence of the 
coherent state in the other mode, we calculate the fidelity F =1 (^f\ Px(t) |*) l5 where 
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Figure 1. Fidelity for the prepared SCS for a = 0.5 (box), 1.0 (cross), 1.5 
(circle), and 2.0 (diamond). Here we used the experimental values T\ = 1 ras 
and T2 = 0.9 ras for the damping time corresponding to mode— 1 and mode— 2, 
respectively. 



pi(t) and |\E r ) 1 is the evolved (mixed) and the prepared (ideal) SCS, respectively, for 
the prepared SCS in mode 1. After a straightforward calculation and rearranging the 
terms, we obtain 

F = A/; 2 AA e 2 e~ (|Q( * )|2+|ao ' 2 ) {cosh [a*a(t) + a a*(t)} 

+ cosh[a* a(t)-a a*{t)]}[2 + Z(t) + Z*(t)], (27) 

where the normalization factors are Af p = [1 + e~ 2 l Q °l 2 ]~5 , W e = [2 + (Z(i)e~ 2 l a OI 2 + 
hc)]~* for the prepared (ideal) and the evolved (mixed) states, respectively, and 
a(t) = a a [un(t) + ui 2 {t)]. 

Fig.l shows the fidelity of the prepared SCS state considering parameters recently 
used in experiments with two modes of a cavity QED [S] . As usual, we have disregarded 
the cavity-field losses during the considerably short time during which the atom crosses 
the cavity, about 10 _4 s. As expected, as the field excitation increases, the fidelity 
decays faster. For a around unity, the fidelity indicates a large time, as compared with 
the damping time, during which the SCS remains available for further operations. 

In conclusion, in this paper we presented, for the first time, a calculation of the 
fidelity for a superposition of coherent states prepared in one of two modes of a cavity 
QED. To achieve this goal, we solved the problem of two modes embedded in a same 
dissipative environment. Different from previous studies, our problem includes off- 
resonantly interaction between a two-level atom, necessary for preparing the coherent 
state superposition, and one of the two modes. Our result indicated high-fidelity of 
the prepared superposition of coherent states when its mean photon number is on 
order of unity. 

We thank the VPG/Vice-Reitoria de Pos-Graduacao e Pesquisa-UCG (NGA), 
CAPES(IPQ, WBC), and CNPq (NGA), Brazilian Financial Agencies, for the partial 
support. 
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